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1 .  INTRODUCTION 


^  In  the  theory  of  electromagnetic  radiation  processes  in  nonequilibrium 
relativistic  beam-plasma  systems,  an  important  role  is  played  by  the  time- 
average  nonlinear  force  acting  on  a  relativistic  test  particle  due  to  the 
total  electromagnetic  field.  A  major  nonlinearity  in  the  field  arises  from 
the  dependence  of  the  Lorentz  force  on  tlie  particle  trajectory,  which,  in 
turn,  depends  on  the  field.  Thereby  the  force  has  effectively  a  nonlinear 
dependence  on  the  total  electromagnetic  field. 

In  the  present  irork,  the  time-average  of  this  nonlinear  force  is  derived 
to  second  order  in  the  total  field  for  the  case  of  unpolarized  and  slightly 
perturbed  relativistic  particle  trajectories.  The  latter  condition  is  that  of 
applicability  of  the  classical  Born  approximation  for  a  plasma,  namely,  that 
the  relativistic  particle  momentum  be  much  greater  than  the  electromagnetic 
impulse  received  by  the  particle  in  a  time  interval  given  by  the  inverse 
plasma  frequency.  Possible  important  screening  effects  arising  from  test 
particle  polarization‘*^re  ignored  here  but  will  be  fully  addressed  in  another 
report.^  The  result  of  the  present  calculation  agrees  with  that  of  Akopyan 
and  Tsytovich."^^  It  is  important  in  calculations  of  collective  radiation  and 
the  conditions  for  the  occurrence  of  radiative  instability  in  beam-plasma 
systems.  ^ _ 

2.  THE  SECOND-ORDER  NONLINEAR  FORCE 

A  charged  test  particle  in  a  beam-plasma  experiences  a  force  given  by 

+  ^dp  * 

Here  a  is  a  label  designating  the  species  of  a  particular  particle,  whose 
position  and  velocity  at  time  t  are  given  by  and  respectively. 

fi(?^,t)  and  6{?^,t)  are  the  electric  and  magnetic  field,  respectively,  at  the 
position  of  the  particle  at  a  particular  time.  The  charge  e^^j  is  the  actual 
charge  of  the  particle.  The  quantity  ^^jp/  the  dynamic  polarization  force,  is 
a  force  correction  arising  from  the  polarization  charge  surrounding  the  test 
particle.  The  latter  arises  from  the  redistribution  of  other  plasma  particles 
due  to  the  induced  fields  around  the  test  particle.  This  polarization  charge 
tends  to  screen  the  actual  charge.  In  the  present  discussion  the  time-average 
of  the  polarization  force  is  ignored;  namely,  it  is  assumed  that  conditions  in 
the  beam-plasma  are  such  that 


<^dp>  •  °  •  (2) 
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It  is  however  important  to  emphasize  that  in  a  complete  analysis  of  collec¬ 
tive  radiation  effects  the  dynamic  polarization  effect  must  be  included.  For 
the  present  discussion,  substituting  equation  (2)  into  equation  (1),  then 
gives  effectively 


The  relativistic  equation  of  motion  governing  the  motion  of  the  particle 
is  given  by 


'a  ^  (^a^a)  =  ' 


vrtiere  m^^^  is  the  rest  mass  of  the  particle  and  is  the  relativistic  gamma 
given  by 


-  [•  -  im 


The  force  may  be  written  in  terms  of  its  Fourier  transform  namely, 


^OL  ® 


-i(u)t-^‘ra) 


where  a>  and  k  are  the  angular  frequency  emd  wave  vector,  respectively,  and 


dk  s  d^kdo)  . 


Similarly,  for  the  electromagnetic  field  Fourier  decomposition,  at  the  parti¬ 
cle,  one  has 

i  .  ,8, 


3  =  /dkS,^  e 


Substituting  equations  (6,8,9)  in  equation  (3),  then 

^odc  =  ®o(®k  • 

Furthermore,  from  the  Maxwell  equation  one  has 
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isees^  wmm  mmwmw  ■awMawa  toumsoa 


It  follows  that 


K  =  oTT-TS-  ' 


where  a  small  imaginary  part  6  of  the  frequency  is  made  explicit, 
ting  equation  (12)  in  equation  (10),  then 


Substitu- 


F  u  =  e 
ok  a 


>  V  X  (k  X 

E,  +  — - r-s - — 

_k  0)  +  i6 


Furthermore,  by  a  well-known  vector  identity. 


5a  «  (If  -  Sfc)  =  (5„-ek)lf  -  (5a-if)«k 


Substituting  equation  (14)  in  equation  (13),  then 


^ok  ~  ®a  ^  y  “  10  +  iSy  ^  to  +  i5 


Then  substituting  equation  (15)  in  equation  (6),  the  force  on  the  particle  to 
second  order  in  the  total  field  is  expressed  in  terms  of  the  particle 
trajectory  and  the  Fourier  component  of  the  total  electric  field  as  follows: 

*  f  ^  -*•  -i(wt-)c»r  ) 

The  time-average  force  acting  on  the  actual  charge  of  the  particle  is  then 
given  by 


<?  >  =  Lim 


ft 


f„(t')df  . 


Substituting  equation  (16)  in  equation  (17)  then 

,  e  r^-/2  /  V  •it  \  ^  \  -i(u)t'-ic»r  ) 

-  ii:  J.t/2 

This  force  depends  on  the  particle  trajectory  through  the  particle  position 
and  velocity  v^  in  equation  (18).  The  classical  Born  approximation  for  a 
plasma  is  to  be  assumed  here,  namely,  that  a  scattering  of  the  particle  by  the 
field  produces  only  a  small  perturbation  of  the  particle  trajectory  in  the 
plasma.  Quantitatively,  one  requires  that  the  particle  momentum  be  much 
greater  than  the  impulse  received  by  the  particle  in  a  time  interval  given  by 
the  inverse  plasma  frequency  (0„,  namely, 

pe 
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Equation  (4)  is  then  to  be  solved  for  the  motion  by  iteration  to  second  order 
in 


Substituting  equation 


-»• 


pe 


I  Pal 


(16)  in  equation  (4)  one  has 

-i  [u)t-ic 


To  lowest  order,  the  motion  of  a  particle  taken  to  be  located  at  the 
coordinates  at  time  t  =  0  and  moving  with  velocity  v^^^  is  given  by 


and 


V  “V 

a  ao 


where  v^^  is  a  constant,  namely,  the  unperturbed  incident  velocity, 
ting  equations  (22,23)  in  the  right-hand  side  of  equation  (21)  and 
ting,  then 


where 


(20) 

(21  ) 

origin  of 

(22) 

(23) 

Substitu- 

integra- 

(24) 

(25) 


and  the  motion  is  unperturbed  in  the  distant  past.  Dividing  both  sides  of 
equation  (24)  by  m„Y^  then 


Using  equation  (5),  one  has  that 


-  *  (l^)]9“' , . 

/J-i^o)  -  +  i6) 


Substituting  equation  (25)  in  the  left-hand  side  of  equation  (33)  and  simpli¬ 
fying,  then 


V  ‘Av 

'oo  oo  I 


^  e  f  -i(“-5t*V(jolt 

~  T^~  ^an*^  - - 

-i(“  -  i^-v^o  + 


Next,  substituting  equation  (34)  in  the  second  term  on  the  right-hand  side  of 
equation  (31),  and  to  the  required  order  replacing  by  y^^,  one  has 


“  LH  +  i6) 

r  \ 

e  V  f  -1  (“-!«•  v  )t 

_  fgXgo.  U  ^  _g - -  . 

Integrating  equation  (35)  with  constants  of  integration  conforming  to  equation 
(22),  and  such  that  the  motion  is  unperturbed  in  the  distant  past,  then  after 
simplifying  one  obtains 


t  =  ^  t  -  - 
a  ao  y 


f  -i  ("-»'•  v  )t 

I  ■*■+  e 

Idk  V  •£ - 

J  oo  K/ 


Y  m  c 
oo  a 


—  V 

2  ao 


(uj  -  k»v  +  i6)^ 

'  rtri  •' 


Equations  (35)  and  (36)  express  the  perturbed  motion  of  the  particle  in  the 
first  Born  approximation.  It  is  this  perturbed  motion  that  enters  in  equation 
(18). 


Substituting  equations  (35)  and  (36)  in  equation  (18),  then 
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where  the  last  exponential  has  been  expanded  to  the  required  order.  Combining 
terras  to  second  order  in  accord  with  equation  (20),  and  noting  that  the  limits 
of  time  integration  are  infinite,  then 


<F  >  =  Lim 


[s,  -  s,  |%j .  it'^_ 


2itie  2  /  '  ^  ^ 

+  ^^‘'*  TV — ;;; —  /dkdk,6(a)  +  <D,-(k  +  k,)»v  ) - - - 

t+«  J  [  (0)  +  i6)(u),  +  i( 


+  +  +  -►  + 

E,  k*v  V  •£. 
k  ao  ao  K. 


(d)  +  i6)(u^  +  i6)((D^  -  +  i^)  (u)  +  i6)(a)^  -  +  i^) 

t  it  \-it,  _ 

(d)  +  i6){dj^  +  i6)  ^  i6)fd)^  +  i6)fd)^  -  +  i6) 

^  ^‘^oto  ^ao'K^ _  K  _ 

c^  (d)  +  i6)fd)^  -  +  16]  (d»  +  i6)fdJi  +  i6)fa)^  -  j^i ‘v^j^  +  i6) 

K  ^ao*®k/“  -  ^*^ao) 


(d)  + 

i6)(d)^ 

t  i6)(d)^ 

- 

i** 

•'T'^ao 

+  i6)2 

iT 

^aoA 

K 

1 

(d)  + 

i6 )  (d)^ 

+  i6)(i»)i 

- 

& 

I  ao 

+  i6) 

oo 

.gk 

(d)  + 

i6 )  (d)^ 

+  i5)(<»i'| 

- 

k  1  •  V 

1  (to 

+  i6)2 

^ujt*V  V  •!  fd)  -  )c*v  1 
k^  ^ao  ao  ao^ 


1?  V  •£.  Jt»v  V  'SL 

(to  K  ao  (to  K  ^ 

”  • 

+  i6)(dj,  -  +  i6)^- 


Next  using  the  property  of  the  (ielta  function  to  replace  by  d)  in  the 

first  integrand,  and  to  replace  (d)^  -  by  -(d)  -  in  the  second 

integrand  of  equation  (40),  and  simplifying,  then  the  latter  becomes 


^•^a^  *^a  +  *^0  ' 


where 


13 


t  u  J  0)  +  i6 


#  (2) 
o 


Lim  e 


/dkdk,  , 

- T-^ - 7  6(<»)  +  0),  - 

(a»  +  i6)  [o)^  +  i6]  ' 


where 


A^J>(k,,k)  - 


\  / 

fs -  fi.  .  +  ^aoi*^i  ~  ^aoi*^1i  _  ^aoi^aoi _ 

^ao”*a  03  -  -  i6  -  16^2  ^  1 

Equation  (43)  cein  be  sy^etrized  k  and  k^  as  follows.  Interchanging  dummy 
integration  variables  (k,u))  with  (k^,(D  ),  and  interchanging  dummy  indices  i 
and  j,  then  equation  (43)  becomes 


#(2)  ^  .  aii  ^ 

«  tiS  ~  « 


or,  equivalently. 


f  (2)  =  Lim  ^  e„ 

“  4.^  t  a 


/■  dkdk,  ,  ,  “" 

J  r;.  i«)K  .  is) *  “I  -  *  '^i)-»<»)ifi'kA,jAS?'(k,k,)  . 


However,  using  equation  (44)  one  has 


6(a)  +  u),  -  (It  +  iti)'v^]A(?’(k,kJ  =  6(u)  +  -  (jt  + 


-  ^r^oo  - 


Using  the  properties  of  the  Kronecker  delta  and  the  Dirac  delta  function 
equation  (47)  may  be  rewritten  as 


6[a)  +  0),  -  (ic  +  it,  +  <*>i  -  (it  +  i^i)*Vao) 


^ooi’^j  -  ^ooj’ 

6.  .  +  - 

L  “  -  lfv_  +  1 


j^li  '^ooi'^ooj 


^ 


("  -  it'Voo  +  i6)2 


Using  equation  (44)  in  equation  (48)  then  the  following  symmetry  relation  is 
obtained: 


6(a)  +  0)^  -  (k  +  it,}*v^^}Aj“Hjt.lt,]  =  «(“  +  -  (k  +  jt,)*v^j^)Af“^*(jt,,k} 

Substituting  equation  (49)  in  equation  (46)  then 


f  dkdk 

iilm  ^ 

t+«>  t  ®a 

J  (0)  +  i6)( 

X  6(u)  +  0)^ 

+ 

1 

Therefore  Fa  can  be  expressed  as  half  the  sum  of  the  two  equivalent  expres¬ 
sions — equations  (43)  and  (50) — namely, 

F^^i  =  Lim  e„  I  +  u  -(it  +  ^,)*Vao)  - - TT  ®ki^  -i 

“  t-H-^  M.  ’  V  aoy  (0)  +  i(5)(a)^+i6}  j5,j 

X  [k  Af“Ult,.lt}  +  it,A!“^*(k^,k}J  . 

Equations  (42)  and  (51)  agree  with  equations  (7)  and  (8)  of  Akopyan  and 
Tsytovich.  ^  It  is  to  be  understood  that  Vq  and  in  that  reference  are 
unperturbed  values  which  are  here  designated  by  Vq^  and  Yg^,  respectively. 
Also,  since  single-wave  particle  resonance  is  ignored  there,  the  principal 
part  of  Af?J*(k,,kJ  in  equation  (51)  is  understood,  in  which  case  j ^  nay  be 
treated  as''  effectively  real  and  becomes  a  common  factor.  The  disparity  of  a 
factor  of  (2ir)“^  and  (2iT)”®  in  equations  (7)  and  (8),  respectively,  of  Akopyan 
and  Tsytovich^  results  merely  from  the  different  Fourier  transform  convention 
chosen  there.  For  example,  in  equation  (5)  of  Akopyan  and  Tsytovich^  the 
Fourier  transform  convention  eiq>loyed  there  has  a  factor  of  (2ir)''^  in  the 
inverse  Fourier  transform  in  the  integration  over  the  three-dimensional  wave 

Ja.  V.  Akopyan  and  V.  N,  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz»  Plazmy,  1  (1975),  673  [Sov,  J.  Plasma  Phys.,  1  (1975),  371}, 


vector  spece,  and  a  factor  of  1  for  the  integration  over  frequency,  giving  a 
total  factor  of  whereas  here  a  total  factor  of  1  is  used  in  equation 

(21).  Therefore,  the  appropriate  factor  for  the  Fourier  transform  itself  is 
(2if)“^  there  and  (2it)“'*  here.  Equations  (41),  (42),  and  (51)  give  the  time- 
average  nonlinear  force  to  second  order  in  the  total  field  for  slightly 
perturbed  relativistic  orbits. 

3.  CONCLUSION 

An  expression — equations  (41),  (42),  and  (51) — has  been  obtained  for  the 
time-average  nonlinear  force  acting  on  an  unpolarized  relativistic  test 
particle  in  a  nonequilibrium  beam-plasma  system.  This  expression  holds  to 
second  order  in  the  total  field  under  the  condition  that  the  classical  Born 
approximation  applies  for  the  plasma.  This  result  has  been  used  in  the  work 
of  Akopyan  and  Tsytovich  in  the  theory  of  collective  bremsstrahlung  in 
nonequilibrium  plasmas. 

The  present  work,  together  with  related  work  by  the  author  '  '  is 
important  for  ongoing  work  in  the  calculation  of  collective  radiation 
processes  and  conditions  for  the  occurrence  of  radiative  instability  in  beam- 
plasma  systems. 
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